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Abstract 

We  consider  an  optimization  problem  in  thermoelectric  cooling.  The 
maximum  achievable  cooling  temperature  in  thermoelectric  cooling  is, 
among  other  things,  affected  by  the  Seebeck  coefficient  profile  of  the 
inhomogeneous  materials.  Mathematically,  the  maximum  cooling  tem¬ 
perature  is  a  non-linear  functional  of  the  Seebeck  coefficient  function. 

In  this  study,  we  solve  this  optimization  problem  exactly. 
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1  Introduction 

In  1822  Thomas  Seebeck  [6]  observed  that  if  two  different  metals  kept  at 
different  temperatures  were  joined,  a  current  would  flow.  In  1834  Jean  Peltier 
[4]  discovered  that  there  is  a  heating  or  cooling  effect  when  electric  current 
passes  through  two  conductors.  It  was  not  until  1851  that  William  Thomson 
(Lord  Kelvin)  [8]  drew  the  connection  between  the  Seebeck  and  Peltier  effects, 
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which  was  the  first  significant  contribution  to  the  understanding  of  thermo¬ 
electric  phenomena.  He  predicted  and  subsequently  observed  experimentally 
that  in  the  presence  of  a  temperature  gradient,  a  single  conductor  with  current 
flow,  will  have  reversible  heating  and  cooling.  With  these  principles  of  ther¬ 
moelectrics  in  mind  and  the  rapid  developments  of  semiconductor  materials 
in  the  late  1950’s,  thermoelectric  cooling  has  become  a  viable  technology  in 
microelectronics  with  applications  in  many  areas  including  flight  vehicles  and 
military  equipments. 

In  thermoelectric  cooling  using  inhomogeneous  materials,  the  maximum 
achievable  cooling  temperature  is  mathematically  given  by  [1]: 


i  S(x)  S(x')dx'dx 

Jj} _ «/0 _ 

o  rL  rx 

/  /  S2(x,)dxl dx 

Jo  Jo 


where  L  is  the  length  of  the  thermoelectric  cooling  element,  and  ZT  is  the 
dimensionless  “figure  of  merit”  [5]  which  puts  a  limit  on  the  maximum  achiev¬ 
able  cooling  temperature  when  a  single  stage  of  uniform  material  is  used.  Since 
1990s,  nanostructured  materials  have  drawn  a  lot  of  attention  because  they  can 
achieve  ZT  values  up  to  2.4  at  room  temperature  [7],  [3],  [9].  In  a  parallel  di¬ 
rection,  a  much  larger  cooling  temperature  beyond  that  of  uniform  materials 
can  be  achieved  by  using  graded  thermoelectric  materials 

In  equation  (1),  x  is  the  length  coordinate  along  the  thermoelectric  cooling 
element,  and  function  S(x)  is  the  Seebeck  coefficient  profile  of  the  inhomoge¬ 
neous  element.  The  Seebeck  coefficient  of  a  material  can  be  varied  by  changing 
the  level  of  doping.  In  semiconductor  production,  doping  refers  to  the  process 
of  intentionally  introducing  impurities  into  an  extremely  pure  semiconductor 
in  order  to  change  its  electrical  properties. 

Without  loss  of  generality,  we  take  L  —  1,  or  equivalently  we  normalize  all 

X 

lengths  by  introducing  x  =  —  and  using  the  dimensionless  x  as  the  independent 
variable.  We  introduce  functional  F[S'(x)]  as 


f  S(x)  f  S{x')dxdx 

qsw]  =  h  ,i ,, -  (2) 

/  /  S2{x')dx'dx 

Jo  Jo 
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In  the  expression  of  F[S'(x)]  defined  in  (2),  the  Seebeck  coefficient  profile  S(x ) 
is  not  allowed  to  be  any  arbitrary  positive  function.  Due  to  the  limitations  in 
manufacturing,  an  acceptable  Seebeck  coefficient  profile  S(x)  must  be  between 
Si  and  S2  (S2  >  Si).  Mathematically,  S(x)  is  restricted  by 

Si  <  S(x)  <  S2,  for  all  x  in  [0, 1]  (4) 

The  goal  of  the  current  study  is  to  find  an  acceptable  Seebeck  coefficient  profile 
that  will  yield  the  largest  maximum  achievable  cooling  temperature.  From 
a  mathematical  point  of  view,  that  is,  to  optimize  functional  F’[S'(x)]  with 
respect  to  function  S(x)  subject  to  constraint  (4). 


2  Exact  solution  of  the  constrained  optimiza¬ 
tion  problem 


To  optimize  functional  F[S'(:r)]  defined  in  (2)  subject  to  constraint  Si  < 
S(x)  <  S2 ,  we  rewrite  the  integral  in  the  numerator  and  the  integral  in  the 
denominator  of  FfS'fT)],  respectively,  as 


S(x)  [  S(x')dx'dx  —  [  [  S{x)S{x')dx' dx 
)  Jo  Jo  Jo 


ax 

ax 


S{x)S{x')dx'  dx  + 
S(x)S(x')dx'dx  + 


r  1  rx 

1 0  Jo 

t  f 

10  Jx' 


S{x)S{x')dx'  dx 
S{x)S{x')dxdx' 


f  [  S{x)S(x')dx' dx 
o  Jo 

1  /  r1  ^  2 
!  S(x)dx 
o 


(5) 


S2(x')dx'dx 


x)dx 


Thus,  the  optimization  problem  becomes 


arg  max 
Si<S(x)<S2 


F[S(x) 


(6) 

(7) 
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where 


F[S(x)' 


1  (/ 

^  f  S2(x)(  1  —  x)dx 


(8) 


In  [2],  based  on  intuitions,  a  Seebeck  coefficient  profile  was  guessed  as  the  solu¬ 
tion  of  optimization  problem  (7).  The  conjectured  optimal  Seebeck  coefficient 
profile  is  given  by  [2] 


Q(x)  =  { 


S\,  0  <  x  <  x\ 

Q  ^  ^ 

- - ,  Xi<X<X2 

1  —  X 

S2,  x2  <  x  <1 


where  q,  x i  and  x2  are  given  by 


Si 

q  =  T 

Xi  =  l- 

x2  =  1  - 


_Q_ 

S  i 

q_ 

s2 


l 

“  2 
=  1  ■ 


Si 

2S2 


(9) 


(10) 


Below,  we  will  show  rigorously  that  the  conjectured  optimal  Seebeck  coefficient 
profile  Q(x)  is  indeed  the  exact  solution  of  the  optimization  problem  (7).  That 
is, 


Q(x) 


arg  max 
Sl<S(x)<S2 


F[S(x) 


(11) 


To  proceed,  we  do  it  in  two  steps: 


•  Step  1:  we  calculate  the  value  of  functional  F[Q(x)]  and  at  the  same 
time  derive  two  properties  of  function  Q(x). 

•  Step  2:  we  use  the  two  properties  derived  in  Step  1  to  prove  that  Q(x)  is 
indeed  the  exact  solution  of  problem  (7).  Mathematically,  we  will  show 

F[S(x)]  <  F[Q(x)]  (12) 

for  all  functions  S(x)  satisfying  Si  <  S(x)  <  S2. 
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Step  1  :  The  integral  in  the  numerator  and  the  integral  in  the  denominator  of 
F[Q(x)]  are  respectively 

r  l  rx\  rx2  q  r  1 

/  Q(x)dx  =  /  Si  dx  +  /  - - dx  +  /  S2dx 

Jo  Jo  Jx  i  1  —  X  Jxo 


xx  1  —  X 
l~X2' 


Sixi  -  gin  (y— +  ‘S’2(1  ~  x 2) 


=  ‘S'i  ( 1  —  -J-  )  —  <?  In 

Ol 

=  S,+<fln('$ 


-fi-f 

<~>2 


V 


1  d 
l~^l)J 


+  s2f 

02 


Si, 

=  9  (  2  +  In  (  ^ 


(13) 


Q2(x)(  1  —  x)dx 


[  Sf(l  —  x)dx  +  [  ( — — — (1  —  x)dx  +  [  ^(l  —  x)dx 

Jo  Jx  1  \1  —  X  J  Jx 2 


T  (1  -  (i  -  x‘)2)  -  ql  111  (ivf )  +  T(1  “ l2) 


1  -  X2' 


S2 


s2 

T 


5X2 


+  g2  In 


'S 


,S  1 


+ 


am 

sfr 

2  s,2 


C  2 

=  y  +  g2ln 
=  g2  f2  + In 


,5! 

S2 

Si 

Substituting  (13)  and  (14)  into  (8)  yields 


F[Q  (x)]  = 


(/  Q(x)dx 


Q2(x)(  1  —  x)dx 
2 

n 

5i, 


rl  2  +  in  ( S 


q2  2  +  In 


'S- 


2  +  ln(^ 


5*1 


Xi 


Multiplying  by  the  denominator  Jq1  Q2(x)(  1  —  x)dx,  we  obtain 


2  + In 


(14) 


(15) 


(I))  l  Q(x)2(l  -  x)dx  -  ^  Q(x)dx^j  =0  (16) 
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Properties  (13)  and  (16)  will  play  a  key  role  in  Step  2  below. 
Step  2:  In  this  step  we  shall  show 


F[S(x)  ]  = 


r-1  N,  2 

S(x)dx 


' o 


S2(x)(  1  —  x)dx 


<  x  (  2  +  In  (  ^ 


(17) 


for  all  functions  S(x)  satisfying  Si  <  S(x)  <  S2.  For  mathematical  conve¬ 
nience,  we  write  S(x)  as  Q(x)  plus  perturbation: 


S(x)  =  Q(x)  +  P(x) 


(18) 


Constraint  Si  <  S(x)  <  S2  on  function  S(x)  implies  the  constraint  below  on 
function  P(x). 


0  <  P(x)  <  S2  —  Si,  0  <  x  <  Xi 
—  (S2  —  Si)  <  P(x)  <0,  X2  <  x  <  1 


(19) 


Note  that  condition  (19)  is  a  consequence  of  condition  (4)  but  (19)  is  not 
equivalent  to  (4).  More  specifically,  (19)  is  weaker  than  (4).  To  prove  (17),  we 
only  need  to  show  that 


G[PM]  =  (2  +  ln(|) 


(Q(x)  +  P(x ))2  (1  —  x)dx 


J  Q(x)  +  P(x)dxsj  >  0 


(20) 

for  all  functions  P(x)  constrained  by  condition  (19).  Expanding  the  squares 
in  (20)  and  using  property  (16),  we  have 

G[P(x)]  =  (2  +  In  ^2  J  Q(x)(l  —  x)P{x)dx  +  J  P2(x)(l  —  x)dxSj 

—2  J  Q(x)dx  J  P(x)dx  —  (/  P(x)dxSj  (21) 

To  further  simplify  G[P(x)],  we  write  Q(x)(  1  —  x)  as 

Q(x)(l  -  x)  =  Q{x)(l  -  x)  -  q  +  q  (22) 

Substituting  (22)  into  (21)  and  using  property  (13)  of  Q(x),  we  arrive  at 

G[P(x)} 
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=  ^2  +  In  ^2  J  (Q(x)(l  —  x)  —  q)  P{x)dx  +  J  P2(x)(l  —  x)dx^j 

—  (/  P(x)dxsj  (23 

It  is  straightforward  to  verify  that  Q(a;)(l  —  x)  —  q  satisfies 


^(l  -  x) 

-  9  >  0, 

0  <  X  <  Xi 

<2(x)(l  -  x)  -  q  =  < 

0, 

Xi  <  X  <  x2 

(24) 

.  52(1  -  2;) 

-  q  <  0, 

x2  <  X  <  1 

Combining  result  (24)  and  constraint  (19)  yields 

>0,  0  <  x  <  x\ 

(Q(x)(l  —  x)  —  q)  P(x)  —  <  0,  x\  <  x  <  x2  (25) 

>  0,  x2  <  x  <  1 

Using  result  (25)  and  the  fact  that  —S2  <  P{x)  <  0  for  x  G  [2:2,1],  we  write 
the  first  term  in  (23)  as 

2  [  (Q(x){l  —  x)  —  q)  P(x)dx  +  [  P2{x){l  —  x)dx 

Jo  Jo 

>  [  (S2(l  —  x)  —  q)  P(x)dx  +  f  P2(x)(l  —  x)dx 

J  X  2  Jo 

—  J  — (1  —  x)j  S2(— P{x))dx  +  J  P2(x){l  —  x)dx 

>  J  —  (1  —  x)^  P2{x)dx  +  J  P2{x){  1  —  x)dx 

=  /  (1  —  x)P2{x)dx  +  [  P2{x)dx  (26) 

Jo  Jx 2  So 


Let  us  introduce  a  new  function: 


R(x)  = 


1  —  x 

S2 

1 

q 


,  0  <  X  <  x2 


x2  <  X  <  1 


We  notice  that  R(x)  is  a  positive  function  in  [0, 1],  and  satisfies 
f1  T~h/  \  7  fX2  1  ,  ,  f1  S2  , 


r  1  rx  2  X  f  * 

/  R(x)dx  =  /  - - da;  + 

Jo  Jo  1  —  2;  Jl-2 
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=  -  ln(l  —  x2)  4  —  (1  -  x2) 

q 

=  l  +  ln(2)  +  ln(j|)  (28) 

Combining  (23)  and  (26),  and  expressing  the  result  in  terms  of  function  R(x), 
we  are  led  to 

™  >  (*  + ■- (|  ))I^-(l!Pix)dX)2  (29) 


To  finish  the  proof,  we  apply  the  Cauchy-Schwartz  inequality  to  (/q1  P{x)dxsj 


< 


Finally,  substituting  (30)  into  (29),  we  conclude 

G[PW]>(1  —  ln(2))/;^>0 


(30) 


(31) 


for  all  functions  P(x)  constrained  by  condition  (19).  It  follows  immediately 
that  function  Q(x)  is  indeed  the  optimal  Seebeck  coefficient  profile  for  maxi¬ 
mizing  the  cooling  temperature. 


3  Conclusions 

In  thermoelectric  cooling,  the  maximum  achievable  cooling  temperature  is 
expressed  as  a  nonlinear  functional  of  the  Seebeck  coefficient  profile  of  the 
inhomogeneous  materials  used.  One  challenge  in  thermoelectric  cooling  appli¬ 
cations  is  to  design  an  optimal  Seebeck  coefficient  profile  so  that  the  cooling 
temperature  is  maximized.  In  manufacturing,  the  Seebeck  coefficient  is  varied 
by  changing  the  level  of  doping  on  a  piece  of  semi-conductor  material.  The 
range  of  the  Seebeck  coefficient  is  limited  so  the  Seebeck  coefficient  profile  is 
constrained  between  two  values.  In  the  study  presented,  we  solved  exactly  this 
constrained  optimization  problem  arised  in  thermoelectric  cooling.  Specifically, 
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we  proved  rigorously  that  a  previously  conjectured  optimal  Seebeck  coefficient 
profile  is  indeed  the  exact  solution  of  the  optimization  problem.  The  methods 
and  techniques  employed  in  the  current  study  may  also  be  useful  for  other 
constrained  functional  optimization  problem. 
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